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Stress Analysis of Multilayer Pressure Vessel 


Shildip D. Urade, Dr.D.V.Bhope, Prof. S.D.Khamankar 


Abstract — Multilayer pressure vessel is designed to work 
under high pressure condition. In this paper, the stress analysis 
of multi-layer pressure vessel made ofa homogeneous and 
isotropic material and subjected to internal pressure is 
considered.Thehoop stresses for 1, 2 & 3-layerpressure vessel is 
calculated theoretically. The modeling of pressure vessel is 
carried out in CATIA V5 and this model is imported in ANSYS 
Workbench where stress analysis is carried out. The shrink fit is 
applied during the CAD modeling of multilayer pressure vessel. 
Both theoretical and FEresults are compared and effect of 
multi-layering on stresses induced and the volume requirement 
to sustain the given pressure is calculated. Also optimization of 
number of layers is carried out for multi-layering the pressure 
vessel. From calculations it is observed that as the numbers of 
layers increases, the hoop stresses decreases. 

Index Terms — Residual Stress, Contact Pressure, 

Multi-layer cylinder, Optimization. 


I. Introduction 

The pressure vessels are used to store fluids under pressure. 
Pressure vessels find wide applications in thermal and nuclear 
power plants, process and chemical industries, in space and 
ocean depths and in water, steam, gas and in air supply 
system, in process Industries like chemical and petroleum 
industries. Design of shrunk-fit multilayer pressure vessel is 
much simplified if each layer is considered as a separate layer 
subjected to the shrinkage pressure and the working pressure. 

To increase the pressure capacity of thick-walled 
layers, two or more layers (multi-layer) areshrunk into each 
other with different diametric differences to form compound 
pressure vessel. In three-layer compound pressure vessel the 
usual practice is to shrink the outer layer 3 on to the 
intermediate layer 2 and then shrink the resulting compound 
pressure vessel on to the inner layer 1 .When the outer layer 
contracts on cooling the inner layer is brought into a state of 
compression.The outer layer will conversely be brought into a 
state of tension. If this compound pressure vessel issubjected 
to internal pressure the resultant hoop stresses will be the 
algebraic sum of those resultingfrom internal pressure and 
those resulting from shrinkage. 


flammable fluids or highly radioactive materials. Because of 
such hazards it is imperative that the design be such that no 

leakage can occur. In addition these vessels have to be 
designed carefully to cope up with the operating temperature 
and pressure.To fulfill these requirements the multilayer 
pressure vessel is designed. For designing the multilayer 
pressure vessel, first analytical analysis is carried out by using 
different mathematical equations. In analytical calculations 
hoop stresses developed in the pressure vessel is carried out. 
Also the total interference and contact pressure is calculated 
at the junction of inner and outer layer. 

To compare and analyze the results the analysis is 
carried out on 1. 2 and 3-layer pressure vessel. The modeling 
of pressure vessel is carried out in CATIA V5 and the same 
model isimported in ANSYS Workbench where stress 
analysis is carried out on CAD model of pressure vessel.To 
optimize the number of layers in multilayer pressure vessel, 
the hoop stress calculation is carried out up to 10-layer 
pressure vessel and the percentage reduction in hoop stresses 
with increase in number of layers is determined. 

For FE and theoretical analysis of pressure vessel, 
hoop stresses are calculated by varying the number of layers. 
For analysis, the internal pressure is assumed as 60 Mpa, inner 
diameter is assumed as 100 mm. The layers of multilayer 
pressure vessel are brought in contact by applying shrink fit 
between the layers ofpressure vessel. The thickness is 
calculated by using the analytical formula and to increase the 
number of layers, thickness is divided equally.Material for all 
the three layers is assumed as same i. e. steel with yield 
strengthof 250 MPa. 


III. Analytical estimation of Hoop stresses in 1 -layer 
Pressure vessel: 

Thickness (t) and outer diameter ( d 2 ) of 1-layer 
pressure vessel is given by equation 1 &2, 


Thickness = t 



1 


(i) 


II. Introduction to Problem, Scope and 
Methodology: 


now, outer diameter (cf 3 ) = d t + (t X 2} (2) 


Pressure vessels often have a combination of high 
pressures together with high temperatures, and in some cases 
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Hoop stress in 1-layer pressure vessel is given by equation 

(3), 


(V-V) 
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IV. Analytical estimation of Hoop stresses in 2-layer Pressure vessel: 

The Hoopstresses induced at outer and inner radius of layer 1 only due to contact pressure(p s i 2 ) at the Junction is given 
by equation (4) and (5), 


ff 01 — 


^e! 2 ( r 2 2 + r l i ) 


r 2 “V 

2^ 


2 ^b 12 ( r 2 ) 


'02 


2 _ r 2 


(4) 


( 5 ) 


The Hoop stresses induced at outer and inner radius of layer 2 only due to contact pressure (p s i 2 ) at the junction is given 
by equations (6) and (7), 


^03 — 


2?b12 C 1 ^) 


tV -r. 


( 6 ) 


_ ^b! 2 (^3 + r 2 Z ) 

^0+ _ r 2 _ r 2 
r 3 r 2 

The Hoop stresses induced due to internal pressure (P;) at inner radius of layer 1 is given by equation (8), 

PiCV + r, 2 ) 


^05 — 


(?) 


tV - r. 


( 8 ) 


The Hoop stresses induced due to internal pressure (Pi) at outer radius of layer 1 and inner radius of layer 2 is given by 
equation (9), 


PLr 


— 


i 1 


r, 2 + r 2 


tV -t’i 


The Hoop stresses induced due to internal pressure (Pi) at outer radius of layer 2 is given by equation (10), 

2-Pi-r, 2 


tj 07 — 


r'q* — Tl 


The resultant hoop stresses at inner radius of layer 1 & 2 is given by the equation (11)& (12), 

P L (r 3 2 + ri 2 ) 2, P sl2 (r 2 z ) 


a 0rl — 


tV - r. 


■V - r. 


(9) 


( 10 ) 


( 11 ) 


°Gr2 — 


Pf r l i 


r, 2 + r, 2 


rv - Th 


Pb 12 ( r 3 2 + r 2 2 ) 


IV - Vn 


( 12 ) 


A. Estimation of Contact pressure (Psn) between layer 1 & 2: 

To obtain optimum values of the contact (shrinkage) pressures P s i 2 which will produce equal hoop (tensile) stresses in both 
thelayers, maximum hoop stresses given by the equations (11) and (12) have been equated. 

Equating equations (11) and (12).i.e. = &Q r 2 and rearranging, 


V + «1 2 

r, 2 /r 3 2 + r 2 2 \l 

~ Pb12 

V+r 2 2 2(r 2 2 ) 1 

r 2 _ r 2 
_ 3 r l 

r 2 2 \ r 3 2 — r l 2 )_ 

, r 3 2 — r 2 Z r 2 Z — r i 2 _ 
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Let, 


Thus, 


Ki 


W-T , 1 



■ sl2 



Ml 

r a a -rt a J 


The total interference between layer 1 & 2 is given by, 


■ r 2 

V + V 

, »2 Z + 'i 2 

E 

„ 2 „ 2 1 

_ 3 r 2 

1 r 2 2 ~ r l 2 . 


(13) 

(14) 


(15) 


(16) 


V. Analytical estimation of Hoop stresses in 3-layer Pressure vessel: 


Figure 1 shows the Radial and Hoop stress distribution in three separate layers of 3-layer pressure vessel. 


r 

3 



Ps£3 



Fig 1. Radial and hoop stress distribution 


Figure 2 and 3 shows the interference between layer 1 & 2 and 2 & 3. Interference is the difference between the outer radius 
of inner layer and inner radius of outer layer.By using this interference the layers are bring in contact with each other. 
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Hoop stresses induced in outer and inner radius of layer 1 only due to contact pressure (P s i 2 ) is given by equations (17) 


& (18), 


a &l ~ P312 




r, ‘ - n 


a 92 ~ 


2-P Sl2 r 2 * 


2 _ ~ 2 


(17) 

(18) 


Radial displacement U rio due to contact pressure at outer radius of layer 1 is given by equation (19), 






(1 -V)( K^X \-V 


2 _ v 2 


(19) 


& ( 21 ), 


Hoop stresses induced in outer and inner radius of layer 2 only due to contact pressure (P s i 2 ) is given by equations (20) 


2P s 12 r 2 2 ^s23 Os" "I” r 2 2 ) 


a03 ~ r 2 _ „ 2 


3 *2 '3 

2 _L ^ 2 


rn ‘ - Tn 


— 


^12 ( r 3 " + r 2 ) ^s23 r 3 ‘ 


rn " —Tn 
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( 20 ) 

( 21 ) 


Radial displacement U r2 idue to contact pressure at inner radius of layer 2 is given by equation (22), 


u r2i = 


r 2 (1 + 7) 


( ft ~ >0 ( !V ?? + v)- (i - to ferrr 


r ? * - r. 


..Tn " - r. 


( 22 ) 


Using equations (19) and (22), total interference 8 i 2 at the contact between layer 1 and 2 is given by equation (23), 

Bu=Url, -Ur,* 


$12 ~ 


r 2 (l-V 2 ') 


fr 3 2 +r 2 2 r 2 2 + ri 2 \ 


323 I « 2 _ „ 2 

V 3 r 2 j 


Radial displacement U r20 due to contact pressure at outer radius of layer 2 is given by equation (24), 

rl-\„ 2 


TJ = — 

E 


2P 3l2 (t - V 2 )r 2 2 l .r a 2 +r 2 2 
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V 
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TV* - Tn 


(23) 


(24) 


(26), 


Hoop stresses induced in outer and inner radius of layer 3 only due to contact pressure (P s i 2 ) is given by equations (25) & 


2 p r 1 
_ zj j 23 f 3 

^05 — 2 „ 2 

r*. — tv 


a G6 ~ 


4 * 3 

^523 (^4^ + r 3 2 ) 

.>1 2 _ 2 

r a * i 


(25) 

(26) 


Radial displacement U r3i due to contact pressure at inner radius of layer 3 is given by equation (27), 


37 


www.erpublication.org 

































International Journal of Engineering and Technical Research (IJETR) 
ISSN: 2321-0869, Volume-2, Issue-9, September 2014 


Vrtt = 


fWaU + l0 




(27) 


Using equations (24) and (27), total interference's at the contact betweenlayer 2 and 3 is given by equation (28), 


&23 ~ U r3i U r2a 


^23 — 


r 3 ( 1 -rq 

L (r. 4 2 + ^ 2 ^ + 

2*Wl 

E 

V r 4 r 3 r 3 ~ r 2 j 

r 3 2 ~r 2 2 _ 


(2S) 


Hoop stress at any radius r in pressure vessel due to internal pressure only is given by, 

2 


®9 = 




2 2 


r 4 

— + 1 

__2 1 
r 


(29) 


After finding hoop stresses at all the radii, the principle of superposition is applied, i.e. the variousstresses are then 
combined algebraically to produce the resultant hoop stresses in the compoundpressure vessel subjected to both shrinkage 
pressures and internal pressure Pi. 


Using equations (29) and (18), maximum hoop stress at the inner surfaces of layer 1 at ri is given by equation (30), 


J 9rl 


= P* 


v+v 


ri — ri 


-2 P 


s 12 


— r-, 


(30) 


Using equations (29) and (21), maximum hoop stress at the inner surfaces of layer 2 at r 2 is given by equation (31), 


p t *i 


J 9r2 


V + r. 


2 _ .y, 2 

t J 


+ 


P&12 2 + r 2 2 ) 2P a23 r 3 


2 2 


(31) 


Using equations (29) and (26), maximum hoop stress at the inner surfaces of layer 3 at r 3 is given by equation (32), 


a 9r3 ~ 




V + r 3 ‘ 


rv‘ - r. 


+ P. 


s23 


V +r 3 ' 


r.“ - 


(32) 


A. Estimation of Contact pressure P sl2 and P s2 3 : 

To obtain optimum values of the contact (shrinkage) pressures P s i 2 and P s23 which will produce equal hoop (tensile) stresses in 
all the three layers, maximum hoop stresses given by the equations (30), (31) and (32) have been equated. 

Equating equations (30) and (31) i.e. a 0r i = a 0r2 and rearranging, 


ffl2 


2r 2 2 . r s 2 + r 2 Z 


r 4 2 + r ± 2 


\r/+r 2 T 

r 2 2 — r± 2 r 3 2 -r 2 2 _ 
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r 4 2 -r^l 


+ P 


2 r 3 


2 _ 
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Where di,d 2 ,d 3 ,d 4 are diameters corresponding to radii r l 5 r 2 r 3 r 4> 
Let, 
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C 3 = 


2r 3 2 


r, * - r. 


Hence equation (33) becomes, 


P a 12 = Pi 


c,\ 


+ p 


s 23 


LC,J 


Equating equation (31) and (32) i. e. a 0r2 = a 0r3 and rearranging, 
„ 2 


■sl2 2 

r 3 


Let, 


c = 


+ r. 


rv * - r. 
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r l 2 (V + r 3 2 ) r l 2 ( r 4 2 + r 2 2 ) 

r 3 2 (r 4 2 -r 1 2 ) r 2 2 (r 4 2 - V) 


„ _ V 
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2r 3 


Hence equation (38) becomes, 


n* -r/ - r, 


^12 = 


C 4 J 




*23 


LQJ 


Equation (37) and (42) have been solved to get P s i 2 and P s2 3 in terms of Pi as follows, 

\C B /C 6 )-(C 2 /C 3 ) 


P 3 12 = Pi 


P&22 ~ Pi 


(c 4 /c 6 )-(c 1 /c 6 ) 

{c 3 /cj- (c 6 /c 4 ) 


(36) 


(37) 


fi 
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(38) 

(39) 

(40) 

(41] 

(42) 


(43) 

(44) 


Putting the values of ti, t 2 and t 3 , the equation (23) and (28) can be written as, 


^12 — 


r 2 (l-V 2 ') 


312 





^23 — 


*3(1-^) 


[p ^ 2+1 

i t22 + 1 V 

2^12 

[ 323 W-i 

\ 2 -l) 

1 

lj 

r-j 

1 

i-i 


(45) 

(46) 


VI. Analytical resultsfor 1 , 2 and 3-layer pressure vessels: 

Analytical results for 1-layer pressure vessel are shown in table no. 1. 

Consider, 

Pj= 60 Mpa, di= 100 mm 

Table 1. Analytical results for 1-layer pressure vessel 
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Parameters 

t (mm) 

d 2 (mm) 

Volume (mm 3 ) 

Hoop stress (ffg) (Mpa) 

Analytical 

14 

128 

5011.44 

247.97 


Analytical results for 2-layer pressure vessel are shown in table no. 2. 
Consider, 


P,= 60 Mpa, d,= 100 mm,d 2 = 114 mm,d 3 = 128 mm. 
Ta ble 2. Analytical results for 2-layer pressure vessel 


Parameters 

P s ±2 (Mpa) 

S 1Z (mm) 

(Mpa) 

<r &z (Mpa) 

Analytical 

2.046 

0.00953 

230.21 

230.21 


Analytical results for 3-layer pressure vessel are shown in table no. 3. 
Consider, 


Pi = 60 Mpa, d|= 100 mm,d 2 = 109.33 mm,d 3 = 118.66 mm,d 4 = 128 mm. 
Ta ble 3.Analytical results for 3-layer pressure vessel __ 


Parameters 

P sl2 

P*2 3 


^Z3 



ff fl3 

Analytical 

1.89 

1.602 

0.00564 

0.0056 

224 

224 

224 


VII. Validation of Analytical results by FEA: 

Using the values of5i 2 &5 23 , shrink fit is applied between layer 1 & 2 and between layer 2 & 3respectively. Contact between 
layer 1 & 2 as well as between layer 2& 3 is applied using contact tool in ANSYS Workbench. The surface between layer 1 & 
2 and 2 & 3 respectively defined as “ROUGH”. 

A. FE analysis of 1 -layer Pressure vessel: 

The data used forCAD modeling is given as follows- 
di = 100 mm, 
d 2 = 128 mm. 



25,00 75.00 


Fig 4. Hoop stresses 


B. FE analysis of 2-layer Pressure vessel: 
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The data used for CAD modeling is given as follows- 
d, = 100 mm, d 2 = 114 mm, d 2 i = 113.981 mm, d 3 = 128 mm. 



C. FE analysis of3-layer Pressure vessel: 

The data used for CAD modeling is given as follows - 

di = 100 mm,d 2 = 109.33 mm, d 2i = 109.319 mm,d 3 = 118.66 mm, d 3i = 118.648 mm, d 4 = 128 mm. 



Fig 7. Contact pressure (Psl2) 


Fig 8. Contact pressure (Ps23) 



Fig 9. Resultant hoop stresses 

VIII. OPTIMIZATION OF NUMBER OF LAYERS IN MULTILAYER PRESSURE VESSEL: 
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From the observations it is concluded that with increase in number of layers, hoop stress reduces. But in detailed study it is 
observed that reduction in hoop stress is not uniform over each layer. 

Thus, it is necessary to know the optimum value of number of layers (n) in designing of pressure vessel to minimize the work 
in multi-layering the pressure vessel, So that the cost of manufacturing of pressure vessel can be reduced. 

Analytical estimation of Hoop stresses in pressure vessel consists of n-layers: 

To calculate stresses of pressure vessel having more than 3-layers some different equation have to use. These equations are 
derived by basic equation to calculate hoop stresses in less time and minimum efforts. 

The equation to calculate hoop stresses in pressure vessel is given as- 

= “ p = + ITT (47) 


The equation given above is to calculate hoop stress due to multi-layering. 
Where, 


r = c 1 xc a xc B . X c r X c r+1 .....x c H 


C r +i is given by. 


J r+1 


Z-g a r4i 

1 +* 3 r+± 


SiK r+1 — 


d r 


( 43 ^ 


IX. Discussion: 

Analytical and FEA (ANSYS) results are summarized in 
Table 4, 5 and 6. 

Table 4. Comparison of analytical and FEA results of 1-layer 
pressure vessel __ 


Parameters 

Hoop stress (ffg) (Mpa) 

Analytical 

247.97 

FEA 

248.14 

Difference (%) 

0.068 % 


Table 5. Comparison of analytical and FEA results of 2-layer 
pressure vessel ___ 


Parameters 

Psl2 (Mpa) 

fffli (Mpa) 

(Mpa) 

Analytical 

2.046 

230.21 

230.21 

FEA 

2.27 

228.26 

231.22 

Difference (%) 

9.8 % 

0.85% 

0.43% 


Table 6. Comparison of analytical and FEA results of 3-layer 
pressure vessel ____ 


Parameters 

Psl2 

(Mpa) 

Ps23(Mpa) 

ff 0i 

(Mpa) 

(Mpa) 

«Tflz(Mpa) 

Analytical 

1.89 

1.602 

224 

224 

224 

FEA 

1.9 

1.7 

224.48 

223.23 

225.78 

Difference 

(%) 

0.52% 

5.76% 

0.21% 

0.34% 

0.78% 


Table 8: Effect of number of layers on hoop stresses 


Sr.no. 

Number of 
layers (n) 

Max. hoop 
stresses (Mpa) 

% Reduction in 
hoop stresses 

1 

1-layer 

247.97 

- 

2 

2-layers 

230.21 

7.16 

3 

3-layers 

224.00 

2.69 

4 

4-layers 

222.019 

0.884 

5 

5-layers 

220.434 

0.71 

6 

6-layers 

219.38 

0.47 

7 

7-layers 

219.00 

0.036 

8 

8-layers 

218.08 

0.57 

9 

9-layers 

217.653 

0.19 

10 

10-layers 

217.40 

0.12 



Number o^ayers 


Fig lO.Effect of Number of layers on hoop stress 

The graph shown in figure 10 represents the effect of 
varying number of layers on hoop stresses. It is showing that, 
as the number of layers increases hoop stresses reduces. 

Also as numbers of layers are increasing the variation in 
rate of change of hoop stress is decreasing. It represents that 
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the effectiveness of layers to reduce the hoop stresses is 
reducing with increasing the number of layers. 

X. Conclusion: 

Analytically, the maximum hoop stress for 1-layer pressure 
vessel is obtained as247.97Mpa, for 2-layer pressure vessel 
230.21Mpa and for 3-layer pressure vessel224 Mpa. 

Thus, after observing these results it is concludedthat - 

• Multi-layering of pressure vessel is very useful in high 
pressure applications. 

• By comparison it is observed that increasing the number of 
layers reduces the hoop stresses in pressure vessel. 

• Multi-layering of vessel decreases the hoop stresses at 
innermost surface of pressure vessel and it can decreases the 
difference between maximum and minimum hoop stress as 
compared to monobloc (1-layer) vessel. 

• In optimization of number of layers (n), after observing all 
the results for hoop stresses on 10 layers it is concluded that 
the multi-layering of pressure vessel is effective up to only 
3-layers. Multi-layering above 3-layers is not advisable as 
percentage reduction in hoop stresses is very less above 
3-layers (n>3). Hence, multi-layering up to 3-layers is 
recommended in designing of multi-layer pressure vessel. 

XI. Nomenclature: 


p 

Internal pressure (Mpa) 

Po 

External pressure (Mpa) 


Hoop stress in the layers (Mpa) 

a 0rl 

Resultant hoop stress in layer 1 (Mpa) 

c 0r2 

Resultant hoop stress in layer 2 (Mpa) 

a 0r3 

Resultant hoop stress in layer 3 (Mpa) 

di 

Inner diameter of layer 1 (mm) 

A 

Outer diameter of layer 1 and Inner diameter of 


layer 2 (mm) 

A 

Outer diameter of layer 2 and Inner diameter of 

a 3 

layer 3 (mm) 

d 4 

Outer diameter of layer 3 (mm) 

r i 

Inner diameter of layer 1 (mm) 

1*2 

Outer radius of layer 1 and Inner radius of layer 2 
(mm) 

r 3 

Outer radius of layer 2 and Inner radius of layer 3 
(mm) 

u 

Outer radius of layer 3 (mm) 

Psi2 

Contact pressure between layer 1 and 2 (Mpa) 

Ps23 

Contact pressure between layer 2 and 3 (Mpa) 

u ri0 

Radial displacement at outer wall of layer 1 (mm) 

u r2i 

Radial displacement at inner wall of layer 2 (mm) 


Total interference at the contact between layer 1 

Ol2 

and 2 (mm) 

u r2 „ 

Radial displacement at outer wall of layer 2 (mm) 

U r3i 

Radial displacement at inner wall of layer 3 (mm) 


Total interference at the contact between layer 2 

023 

and 3 (mm) 

V 

Poisson’s ratio 

K r+ i 

Ratio of outer diameter to inner diameter of 
(r+1 )* shell (dr+1/dr) 
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